Brazil-nut effect versus reverse Brazil-nut effect in a moderately dense granular fluid 
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A new segregation criterion based on the inelastic Enskog kinetic equation is derived to show the 
transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE) by varying 
the different parameters of the system. In contrast to previous theoretical attempts the approach is 
not limited to the near-elastic case, takes into account the influence of both thermal gradients and 
gravity and applies for moderate densities. The form of the phase-diagrams for the BNE/RBNE 
transition depends sensitively on the value of gravity relative to the thermal gradient, so that it is 
possible to switch between both states for given values of the mass and size ratios, the coefficients 
of restitution and the solid volume fraction. In particular, the influence of collisional dissipation on 
segregation becomes more important when the thermal gradient dominates over gravity than in the 
opposite limit. The present analysis extends previous results derived in the dilute limit case and is 
consistent with the findings of some recent experimental results. 
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Segregation and mixing of dissimilar grains in agitated 
granular binary mixtures is one of the most important 
problems in granular matter both from a fundamental 
and a practical point of view. In some cases it is a desired 
and useful effect to separate particles of different types, 
while in other processes it is undesired and can be diffi- 
cult to control. Usually, the larger intruder particles tend 
to climb to the top of the sample against gravity (Brazil- 
nut effect, BNE), but under certain conditions they can 
also accumulate at the bottom (reverse Brazil-nut effect, 
RBNE). However, although there is an extensive observa- 
tional evidence of these phenomena, much less is known 
on the physical mechanisms involved in this problem [l[ . 
Among the different competing mechanisms proposed to 
explain the transition BNE<^> RBNE @, H, 0, H , thermal 
(Soret) diffusion becomes the most relevant one at large 
shaking amplitude where the system resembles a granu- 
lar gas and kinetic theory becomes useful to study seg- 
regation. Some previous theoretical attempts have been 
reported in the literature analyzing thermal diffusion ef- 
fects on segregation in bi-disperse granular systems. Nev- 
ertheless, these early contributions have been restricted 
to elastic @ and quasielastic particles 0, have consid- 
ered thermalized systems (and so the segregation dynam- 
ics of intruders is only driven by the gravitational force) 
[U 0], and/or have been limited to dilute gases [1, EJ 
The main goal of this paper is to propose a theory based 
on a recent solution of the inelastic Enskog equation [ll| 
that covers some of the aspects not accounted for in pre- 
vious works. The theory subsumes all previous analyses 
d, 0, B @| j which are recovered in the appropriate lim- 
its. Furthermore, the theoretical predictions are in qual- 
itative agreement with some molecular dynamics (MD) 
simulation results 0, d, [13] and are also consistent with 
previous experimental works [1, [B| . 

We consider a binary mixture of inelastic hard disks 
(d = 2) or spheres (d = 3) where one of the components 
(of mass toq and diameter oq) is present in tracer con- 
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FIG. 1: (Color online) Sketch of the problem studied here. 
The small circles represent the particles of the granular fluid 
while the large circles are the intruders. The BNE (RBNE) 
effect corresponds to the situation in which the intruder rises 
(falls) to the top (bottom) plate. 



centration. In this limit case, one can assume that (i) 
the state of the dense gas (excess component of mass m 
and diameter a < ao) is not affected by the presence of 
tracer particles and, (ii) one can also neglect collisions 
among tracer particles in their kinetic equation. This 
is formally equivalent to study an intruder in a dense 
granular gas, and this will be the terminology used here. 
Collisions among gas-gas and intruder-gas particles are 
inelastic and are characterized by two independent (con- 
stant) coefficients of normal restitution a and ao, respec- 
tively. The system is in presence of the gravitational field 
g = —ge z , where g is a positive constant and e z is the 
unit vector in the positive direction of the z axis. To flu- 
idize the mixture and reach a steady state, in most of the 
experiments energy is injected into the system through 
vibrating horizontal walls. Here, in order to avoid the use 
of vibrating boundary conditions, particles are assumed 
to be heated by a stochastic-driving force which mim- 
ics a thermal bath. Although previous experiments fl4j 
have shown a less significant dependence of the temper- 
ature ratio Tq/T on inelasticity than the one obtained in 
driven steady states fl3j |. the results derived in Ref. 
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FIG. 2: (Color online) Phase diagram for BNE/RBNE for = 
0.1 in the absence of gravity for two values of the (common) 
coefficient of restitution a. Points above the curve correspond 
to A > (BNE) while points below the curve correspond to 
A < (RBNE). The dashed line is the result obtained for 
a — 0.8 assuming energy equipartition (To = T). 



for To/T from this stochastic driving method compare 
quite well with MD simulations of agitated mixtures [B[ . 
This agreement suggests that this driving method can be 
seen as a plausible approximation for comparison with 
experimental results. A sketch of the geometry of the 
problem is given in Fig. [T] 

The thermal (Soret) diffusion factor A is defined at 
the steady state with zero flow velocity and gradients 
only along the vertical direction [z axis). Under these 
conditions, the factor A is defined by [j| 



Ad, InT = d r In 



(i) 



where uq and n are the number densities of the in- 
truder and the fluid particles, respectively. Let us assume 
that gravity and thermal gradient point in parallel direc- 
tions (i.e., the bottom plate is hotter than the top plate, 
d z InT < 0) (see Fig. [TJ. Obviously, when A > 0, the in- 
truder rises with respect to the fluid particles (BNE, i.e., 
d z ln(n /n) > 0) while if A < 0, the intruder falls with re- 
spect to the fluid particles (RBNE, i.e., d z \n(n /n) < 0). 

In order to determine the dependence of the coefficient 
A on the parameters of the system, we focus our attention 
on an inhomogeneous steady state with gradients along 
the z direction. Since the flow velocity vanishes, the mass 
flux of the intruder j z = 0. In addition, the momentum 
balance equation yields 

where p is the pressure and p — ran is the mass density of 
the fluid particles. Upon writing Eq. ([2|), we have taken 
into account tha t p depends on z through its dependence 
on n and T 



11] 



To first order in the spatial gradients 
(Navier-Stokes description), the constitutive equation for 
the mass flux of the intruder is 



11] 
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FIG. 3: (Color online) Phase diagram for BNE/RBNE for 
a = 0.8, \g*\ — 1 and three different values of the solid volume 
fraction d>. 




FIG. 4: (Color online) Phase diagram for BNE/RBNE for 
cj> = 0.5 in the absence of thermal gradient (\g*\ — * oo) for 
two values of a. The dashed and dashed-dotted lines refer to 
the results obtained by Jenkins and Yoon [6| and Trujillo et 
al. [7|, respectively. 



where Do, D, and D T are the relevant transport coeffi- 
cients. Expressions for the pressure p and the transport 
coefficients Dq, D, and D T have been recently obtained 
in the undriven case by solving the Enskog kinetic equa- 
tion by means of the Chapman-Enskog method in the 
first Sonine approximation (llT | . The extension of these 
results to the driven case is straightforward. The con- 
dition j z — along with the balance equation ([2|) allow 
one to get the thermal diffusion factor A in terms of the 
parameters of the mixture as 



A = 



(3D T * - (p* + g*){D* a + D*) 



(4) 



Here, f3 = p* + tfid^p*, g* — pg/nd z T < is a dimen- 
sionless parameter measuring the gravity relative to the 
thermal gradient, and p* = p/nT = 1 + 2 d_2 x0(l + ot). 
Moreover, x('/ > ) is the pair correlation function for the 
granular gas, </> = [n d/2 /2 d - 1 dT(d/2)]na d is the solid vol- 
ume fraction and the reduced transport coefficients are 
explicitly given by 
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Here, 7 = Tq/T is the temperature ratio, M = mo/m is 
the mass ratio, cu = erg /cr is the size ratio and 



^-D = ^77^ X01 



' 7 + M 
M(l + M) 



(l + oo). (8) 



In addition, xo is the intruder-gas pair correlation func- 
tion and /Lto is the chemical potential of the intruder. 
When the granular gas is driven by means of a stochas- 
tic thermostat, the temperature ratio 7 is determined [9] 
from the requirement 7^0 = MQ where the expressions 
for the cooling rates £0 an d C m the local equilibrium 
approximation can be found in Ref. [111 ]. 

The condition A = provides the segregation cri- 
terion for the transition BNE O RBNE. Since (3 and 
Dq are both positive, then according to ((4]), sgn(A) = 
sgn(/3D T * - (p* + g*){D* a +£>*)). Consequently, taking 
into account Eqs. ©-((T]), the segregation criterion can 
be written as 

.9*(7 " Mff) - 7^ + M ^ Xo0(l + <*o) 



M + -v 



2 1 + AT 
= 0, 



(9) 



where A = [(1 + Lu)- d /T X o}(d<plio)T,n - Equation © 
gives the phase-diagram for the BNE/RBNE transition 
due to Soret diffusion of an intruder in a moderately 
dense granular gas. This is the most relevant result of 
this paper. The parameter space of the problem is six- 
fold: the dimensionless gravity g*, the mass ratio mo/m, 
the ratio of diameters <tq/<j, the solid volume fraction </>, 
and the coefficients of restitution a and a^. The influ- 
ence of density on segregation is accounted for by the 
second and third terms in Eq. ([9]). As expected, when 
niQ = m, <7q = a and a = ctQ, the system (intruder plus 
gas) is monodisperse and the two species do not segre- 
gate. This is consistent with Eq. © since in this limit 



case D 1 * = Dn + D* = and so, the condition 



for any value of a and 4 
(3 = p* = 1 so that Eq. 



In the dilute limit case 
M reads 



holds 

-0), 

(10) 



«?*( 7 -M) = 0. 



Note that in the absence of gravity (g* = 0), Eq. p0|) 
applies for any value of the parameters of the system 
and so the intruder does not segregate in a dilute gas. 
When g* ^ 0, the solution to JTO]) is 7 = M. This result 
agrees with the recent segregation condition derived from 
the Boltzmann equation 



In general, segregation is driven and sustained by both 
gravity and temperature gradients. The combined effect 
of both parameters on Soret diffusion is through the re- 
duced parameter g*. This parameter measures the com- 
petition between these two mechanisms on segregation. 
In previous theoretical studies on dense gases (g, l7(, the 
temperature was assumed to be uniform in the bulk of 
the system (d z T = 0) so that segregation of intruder was 
essentially driven by gravity. This is quite an interest- 
ing limit since most of the experiments d H GJ] have 
been carried out under these conditions. In this limit 
— > 00), Eq. ([9]) becomes 
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while the segregation criterion found independently by 
Jenkins and Yoon [U (for an elastic system) and by Tru- 
jillo et al. 0] is 
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(12) 



Equation (fTTjl reduces to Eq. (fT2|) when one (i) neglects 
the dependence on inelasticity and assumes equipartition 
in certain terms, (ii) takes the approximation A = 1 
(which only applies for a dilute gas of mechanically equiv- 
alent particles), and (iii) neglects high density corrections 
(last term in the denominator of (|11|) ). Thus, even in the 
particular limit \g*\ — > 00, the criterion (jlll) is much 
more general than the one previously derived [7[ since it 
covers the complete range of the parameter space of the 
problem. 

Henceforth and for the sake of concreteness, we assume 
a three-dimensional system with a — ao- In this case 
[Rl,X = (l-^)/(l-0) 3 and 
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The expression for the chemical potential consistent with 
the above approximations can be found in Ref. [T3|. A 
phase diagram delineating the regimes between BNE and 
RBNE is shown in Fig. [2] for <f> = 0.1, g* = and two 
values of the (common) coefficient of restitution a. We 
observe that, in the absence of gravity, the main effect 
of dissipation is to reduce the size of the BNE. It is also 
apparent that the RBNE is dominant at both small mass 
ratio and/or large diameter ratio. In addition, compari- 
son with the results obtained for a = 0.8 assuming that 
Tq = T shows that the nonequipartition of granular en- 
ergy has an important influence on segregation in the 
absence of gravity. This is consistent with MD-based 
findings of Galvin et al. [l2j. Moreover, the results ob- 
tained from Eq. © show that the form of the phase di- 
agrams depend significantly on the value of the reduced 
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gravity \g*\ (namely, reverse buoyancy relative to the ef- 
fect of Soret diffusion). Thus, in general for given values 
of m /m, do/ a, ag, a and <f>, a transition between both 
states BNE/RBNE is possible by changing \g*\. To il- 
lustrate the influence of the reduced gravity, a phase dia- 
gram is plotted in Fig. [3] when \g*\ — 1 (gravity compara- 
ble to the thermal gradient) for different values of the vol- 
ume fraction c/>. In contrast to Fig. [21 it is apparent that 
the RBNE regime appears essentially now for both large 
mass ratio and/or small size ratio. With respect to the 
dependence on density, Fig. [3] shows that in general the 
role played by density is quite important since the regime 
of RBNE decreases significantly with increasing density. 
Following Trujillo et al. 7], since the effect of shaking 
strength of vibration on the phase diagram BNE/RBNE 
can be tied to the effect of varying the volume fraction <f>, 
it is apparent from Fig. [3] that the possibility of RBNE 
will increase with increasing shaking strength. This fea- 
ture agrees with the experimental findings of Brcu et al. 
[3| since their results show similar behavior with the ex- 
ternal excitation. The form of the phase diagram in the 
limit |<7*| — > oo is shown in Fig. U for (f> = 0.5 and two 
values of the coefficient of restitution a = 1 and 0.9. Our 
results indicate that, in contrast to the case of Fig. [21 the 
main effect of inelasticity is to reduce the size of RBNE 
region, which qualitatively agrees again with experiments 
On the other hand, our predictions disagree with the 
theoretical results derived by Trujillo et al. Q since the 
latter found that the mass ratio is a two-valued function 
of the size ratio and so the main effect of dissipation is to 
introduce a threshold value of the size ratio above which 
there is no RBNE. The results also indicate (not shown 
in Fig. 3]) that nonequipartition has a weaker influence 
on segregation when \g*\ — > oo than in the opposite limit 
(1.9*1 = 0). This behavior qualitatively agrees with the 
experiments carried out by Schroter et al. [5j for vibrated 
mixtures since they find energy nonequipartition to have 
no discernible influence on their results. 

In summary, a kinetic theory based on a solution of 
the inelastic Enskog equation has been used to analyze 
thermal (Soret) diffusion effects on segregation for an in- 
truder in a driven moderately dense granular gas under 
gravity. The present study goes beyond the weak dissipa- 
tion limit, takes into account the influence of both ther- 
mal gradient and gravity on thermal diffusion and applies 
for moderate densities. Although the theory is consistent 
with previous numerical and experimental results, a more 
quantitative comparison with the latter would be desir- 
able. As a first test, kinetic theory predictions in the 
Boltzmann limit Q compare well with MD simulations 
of agitated dilute mixtures [8(. Given that the results 
derived here extends the description made in Ref. [9J to 
moderate densities, it can be reasonably expected that 
such a good agreement is also kept at finite densities. In 
this context, it is hoped that this paper stimulates the 
performance of such simulations. On the other hand, it 



must be stressed that the present work only deals with 
the tracer or intruder limit. This precludes the possibil- 
ity of comparing our theory with the results reported by 
Schroter et al. [eJ in agitated mixtures constituted by 
particles of the same density and equal total volumes of 
large and small particles. When convection is practically 
suppressed, they studied the influence of dissipation on 
Soret diffusion. I plan to extend the present theory to 
finite concentration to carry out a comparison with the 
above computer simulation results H in the near future. 
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